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In this paper we study information flow paths in a data network, where traffic 
generated by servers ( or sources) takes a multi-hop path in order to reach its 
clients (destinations). Each node in the middle of this multi-hop path should 
route the incoming traffic and the traffic generated by itself to the next hop 
in such a way that the traffic reaches its destination while avoiding congestion 
in the links. For simplicity, we will assume the network only carries single 
commodity traffic, i.e., all of the traffic should be routed to a single destination. 

Let us assume that the data network is modeled by a connected graph 
G(V,E) where the set of vertices, V (\V\ = N) represents the nodes (or data 
servers) and the set of edges, E (\E\ = M) represents the links. For our analysis, 
we assume that the destination for all of the traffic in the network is node N. 
Two nodes are called adjacent if they are connected by an edge. A cut-set is a 
set of edges that if they are removed from the graph it becomes disconnected. 
A minimum cut-set is a cut-set with minimum cardinality. The throughput of 
a network with arbitrary sources and destinations is limited by the capacity of 
the links in the minimum cut-set of the network. If traffic is not routed care- 
fully, when the network is under heavy traffic, some of these links may become 
congested while other are under-utilized. It is desirable to route the traffic such 
that each link in the minimum cut-set achieves its maximum capacity. When the 
capacity of all of these links are the same, an optimal routing would route the 
traffic on the bottleneck links such that the traffic is distributed evenly among 
them. 

We arbitrarily assign a direction to each edge e m such that if it connects 
nodes r and s where r < s, then the the direction of e m is from r to s. Let 
us denote the traffic load on edge e m with I m . Note that if traffic flows from 
node r to node s then I m is positive and it is negative if traffic flows in opposite 
direction. For each node n, the sum of incoming traffic load (from adjacent 
nodes with lower index ) and the traffic generated by the node itself (T n ) is 
equal to the sum of outgoing traffic (to adjacent nodes with higher index). 
These conservation of traffic constraints are written for nodes 1 ... N — 1 in 
matrix form as 

C-I = T (I) 

where C is a (N — 1) x M matrix with elements in {—I, 0, +1 }, I is a M x 1 
vector of traffic loads I m and T is a iV - 1 x 1 vector of traffic sources T n . For 
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— 1 if ej connects node i to a node with lower index 
lif ej does not start or end with node i (2) 
+ 1 if ej connects node i to a node with higher index 

We define a cost function 

M 

mi) = \m = e p m ( 3 ) 

and write the following optimization problem: 

minimize Jp(I) 

(4) 

subject to C ■ I = T 

In order to solve this problem in a distributed fashion we employ Sequential 
Quadratic Programming (SQP) in conjunction Jacobi iterations. 

Let i be a small perturbation in I such that C ■ (I + i) = T. Thus C ■ i = 
. We also have 

M 

771=1 

^J p {I) + -i T -Q-i + S T -i (6) 

where S = p^' 1 , P' 1 ) and Q = diag (p(p - 1)(7?- 2 , . . . j£r 2 )) . Thus, 
for a given initial value of I, problem (JTTJ becomes 

minimize </»(«) = -« T ■ Q ■ i + S T ■ i 

I pw 2 ^ (7) 

subject to (7 • i = 
We can write Jp(i) as 

M 

M*) = E = 2 P(P - ^^m+Pl^m (8) 
m— 1 

= E>-w- 2 (« m+ ^) 2 -(^ p ) do) 

defining a new variable i m = i m + ^fy we have the following optimization 
problem 
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minimize J p (i) = r m i 2 m 

1 m=l (11) 

subject to C • i = — !— CT 

P- 1 

where r m = \p{p — 1)I^~ 2 . The problem in this canonical form is easy to solve 
distributedly using Jacobi iterations. For each node we define a potential u such 
that if nodes r and s with r < s are connected by edge m then we have 

lm = (12) 

I'm 

Each node calculates its potential only as a function of the potentials of its 
adjacent nodes. Let A n be the set of nodes that are adjacent to node n. Let 
£ n (fc) be the index of the link that connects node n to an adjacent node k. Then 
by Jacobi iterations we have 



p-i Z^ keA rtn(fc) (13) 



V" 1 

. Note that the potential for the destination (node N) is always zero. After 
calculating the potentials for all of the nodes it is i m is 



p-1 



(14) 



To summarize, up to here we calculated distributedly a correction vector i 
so that changing the traffic loads from I to I + i reduces the cost function J P (I). 
Since J p is a convex function of I for p > 1, this SQP steps can be repeated 
several times in order to reach the global minimum of J p . Note that as p — > oo 
the optimization problem becomes a minimax problem trying to distribute the 
traffic load evenly in the minimum cut-set links. 
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